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Abstract 

The e + e~— > e + e~hadrons reaction, where one of the two electrons is detected 
in a low polar-angle calorimeter, is analysed in order to measure the hadronic pho- 
ton structure function F% . The full high-energy and high-luminosity data set, col- 
lected with the L3 detector at centre-of-mass energies 189 GeV < ^/s < 209 GeV, 
corresponding to an integrated luminosity of 608 pb _1 is used. The Q 2 range 
11 GeV 2 < Q 2 < 34 GeV 2 and the x range 0.006 < x < 0.556 are considered. 
The data are compared with recent parton density functions. 
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1 Introduction 



Photons are ideal tools for probing the structure of more complex objects such as the proton 
in deep- inelastic scattering experiments. At LEP, in the e + e~ — > e + e~7*7 — > e + e~ hadrons 
reaction, two virtual photons are produced by the incoming electrons 1 ) and their interaction 
yields hadrons. If the scattering angle of one of the electrons, # t ag, is sufficiently large, it 
is observed in the low polar-angle electromagnetic BGO calorimeter [1] of the L3 detector [2], 
originally devised to detect low angle Bhabha scattering in order to measure the LEP luminosity. 
This allows to measure the four-momentum, k' ', of this "tagged" electron. For "single-tagged" 
events the second electron is undetected, its polar angle is small and the virtual photon radiated 
from this electron is quasi-real. In the framework of a deep-inelastic scattering formalism the 
process e + e~— > e + e~ hadrons is written as the convolution of the target photon flux with the 
reaction e(k) + ^(p) — > e(k') + hadrons. The photon, 7*, with four-momentum q = k — k' 
and a large virtuality Q 2 = —q 2 ~ 2E tllg E beilTa ( y l — cos6> tag ), is considered as a point-like probe 
investigating the structure of the target photon, 7, with four-momentum p and virtuality P 2 = 
—p 2 ~ 0. E tag is the energy of the tagged electron and -Ebeam the energy of the beam. The 
differential cross section is written in terms of the scaling variables x = Q 2 /2(p- q) = Q 2 /{Q 2 + 
W 2 ^ + P 2 ) and y = (q- p)/(k ■ p) = 1 - (£ tag /£ bcam cos 2 tag ) as [3,4]: 

dcr e ^^ e x (x, Q ) _ 27TQ ^ | ^ „.\2\jry/„ /^2\ „ 2 7717/ s^2\ 



dxdQ 2 xQ l 



-[(1 + (1 - y) 2 )F2(x, Q 2 ) - y 2 F2(x, Q 2 )] (1) 



The variable x depends on the two-photon centre-of-mass energy, VF 77 , equal to the effective 
mass of the produced hadrons. The inelasticity y is small (y < 0.3) in the kinematic region of 
this study and consequently only F^x^Q 2 ) contributes appreciably to the cross section. By 
convention, F^jot is measured, where a is the fine-structure constant. Using this approach, the 
photon structure function has been extensively studied at low-energy e + e~ colliders [5] and at 
LEP [5-7]. 

A virtual photon can interact as a point-like particle in "direct processes" ; it can fluctuate 
into a vector meson (p,u,<f>) inducing soft hadronic interactions in "VDM processes" or it 
can interact via its partonic content of quarks or gluons in "resolved processes". High Q 2 
single-tag events favour perturbative QED and QCD diagrams such as 77 — > qq , 7q — > gq or 
7g — > qq. The two resolved processes, 7q — > gq and 7g — > qq, are described using parton density 
functions extracted from the photon structure functions measured in previous experiments at 
PEP, PETRA and TRISTAN. Reviews of the existing parameterisations may be found in 
References 8 and 9. Recently, a new parametrisation was obtained adding published LEP 
data [10]. 

This analysis uses the 608.1 pb -1 of high-energy LEP data, collected at e + e~ centre-of- 
mass energies 189 GeV < y/s < 209 GeV. The data are grouped in four average y/s values, 
presented in Table 1. These high-energy data allow to extend our previous measurements [6] 
at ^/s ~ 91 GeV and ^/s = 183 GeV in the small- a; region down to 0.006 and in the medium-a; 
region up to 0.556 for the Q 2 range 11 GeV 2 < Q 2 < 34 GeV 2 . The e + e~— > e + e~hadrons cross 
section and the K7 photon structure function are studied as a function of x in the three Q 2 
intervals 11 GeV 2 < Q 2 < 14 GeV 2 , 14 GeV 2 < Q 2 < 20 GeV 2 and 20 GeV 2 < Q 2 < 34 GeV 2 . 
The Q 2 evolution of F 2 7 is also studied combining the values at < Q 2 >= 12.4, 16.7 and 
25.5 GeV 2 with our previous measurements. 



Throughout this Letter, the term "electron" indicates both electron and positron. 
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2 Monte Carlo Models 



The value of the Q 2 variable is accurately determined by measuring the four-momentum of 
the scattered electron. However, the effective mass of the final state hadrons is only partially 
reconstructed, as these are often produced at low polar angles where no tracking system can 
be installed. A Monte Carlo modelling of the final state hadrons is therefore necessary [11] to 
determine the x variable. 

Three Monte Carlo generators are used to model the process e + e~— > e + e~ hadrons: PHO- 
JET [12], PYTHIA [13] and TWOGAM [14]. 

PHO JET describes hadron-hadron, photon-hadron and photon-photon collisions. It is based 
on the Dual Parton Model combined with the QCD-improved parton model [15]. In order to 
have a continuous transition between hard and soft processes, the distribution of the transverse 
momentum, p t , of the soft partons is matched to the one predicted by QCD. The two-photon 
luminosity is calculated from the flux of transversely polarised photons; corrections for the lon- 
gitudinally polarised photons are thus incorporated into an effective two-photon cross section. 
The transition from real-photon to virtual-photon scattering is obtained by a change of the 
relative weight of all partial cross sections. 

PYTHIA is a general purpose Monte Carlo. For two-photon interactions it incorporates 
leading order (LO) hard-scattering processes as well as elastic, diffractive and low p t events. 
The classification of the photon interactions into three different components, direct, resolved 
and VDM, results in six different classes of events. Events are also classified according to the 
hard scales involved in the process: photon virtualities and parton transverse momenta. 

TWOGAM generates three different processes separately: point-like photon-photon inter- 
actions, resolved processes, and non-perturbative soft processes described by the Generalised 
Vector Dominance Model (GVDM). The structure of the program is modular and the photon 
flux is calculated with an exact LO formula. The cross sections of the three different processes 
are adjusted to fit the ^-distribution of the data. The cross section of the direct process is 
fixed to the value expected in our kinematic range, a = 41 pb. The QCD and the VDM cross 
sections are then adjusted to a = 5 pb and a = 28 pb, respectively. 

For the three Monte Carlo generators parton showering and hadronisation are described by 
JETSET [16]. The dominant backgrounds are evaluated with PYTHIA for e + e~ — > qq(7) and 
DIAG36 [17] for e+e" -> e+e'r+r-. 

All Monte Carlo samples are generated with a luminosity at least five times greater than 
the experimental one. All events are passed through a full detector simulation which uses the 
GEANT [18] and GHEISHA [19] programs and takes into account detector efficiencies and 
time-dependent effects. Monte Carlo events are then reconstructed in the same way as the 
data. 

3 Data Analysis 

Events are mainly accepted by two independent triggers: the single-tag trigger and the central 
track-trigger. The single-tag trigger requires at least 70% of the beam energy to be deposited 
in one of the low polar-angle calorimeters, in coincidence with at least one track in the central 
tracking chamber. The central track-trigger requires at least two tracks back-to-back in the 
transverse plane within ±60°, each with p t > 150 MeV. The average trigger efficiency is about 
97%. There is a single-tag trigger signal for about 90% of the selected events and a central 
track-trigger signal for about 85% of the selected events. 
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3.1 Event Selection 



Events are selected by requiring a single scattered electron in the low polar- angle calorimeter 
and a hadronic final state. A tagged electron candidate is the highest energy cluster with a 
shape consistent with an electromagnetic shower, E tag / E^ eSum > 0.7, as shown in Figure la, 
and a polar angle in the fiducial region 0.0325 rad < 9 < 0.0637 rad inside the geometrical 
acceptance 0.030 rad < 9 < 0.066 rad. To ensure that the virtuality of the target photon 
is small, the highest-energy cluster in the low polar-angle calorimeter opposite to the tagged 
electron must have an energy less than 20% of the beam energy, as shown in Figure lb. 

At least four additional particles must be detected. A particle can be a track or a photon. 
A track must have p t > 100 MeV and a distance of closest approach in the transverse plane to 
the interaction vertex of less than 10 mm. A photon is a cluster in the electromagnetic BGO 
calorimeters with energy above 100 MeV, not associated with a charged track. 

To reduce the background from the process e + e~ — > qq(7), the total energy deposited in 
the electromagnetic and hadronic calorimeters must be less than 40% of the center-of-mass 
energy, as shown in Figure lc. The events with a large value of the total energy are due to the 
e + e~ — > Z7 — > qq7 process, where the radiative photon is misidentified as the tagged electron. 

The mass of the hadronic final state, W vis , is calculated from all tracks and calorimetric 
clusters. Additional clusters detected in the low polar-angle calorimeter are assigned the pion 
mass and included in the calculation of W V1S . To avoid the hadronic resonance region, W V - 1S is 
required to be greater than 4 GeV, as presented in Figure Id. 

Figure 2 shows the Q 2 distribution for each y/s sample. Only events with 11 GeV 2 < 
Q 2 < 34 GeV 2 are studied. The number of selected events and the backgrounds from the 
e + e~ — > e + e~T + r" and e + e~ — > qq(7) processes are given in Table 1. The background is 
dominated by the e + e~ — > e + e~r + r _ production. The contribution from the e + e~ — > t + t~ 
and e + e~ — > W + W~ processes is negligible. The background from beam-gas and beam-wall 
events is found to be negligible by inspection of the radial distribution of track intersections. 
The W vis and x vis = Q 2 / (Q 2 + W 2 is ) distributions are presented in Figure 3 for all selected data. 
The PYTHIA and TWOGAM model reproduce the data rather well, except at large values of 
W vis . PHOJET presents a harder mass spectrum and predicts too many events for x vis < 0.1 
and is therefore not used in the following. 

The total acceptance is calculated for each data sample separately. It takes into account the 
trigger efficiency, the geometrical acceptance and the selection cuts. An example is presented 
in Figure 4a for the data at y/s— 189 GeV. 

4 Results 

4.1 Unfolding and Differential Cross Sections 

The cross section Acr/Ax as a function of x for the reaction e + e~— > e + e~hadrons is measured 
for three Q 2 intervals: 11 GeV 2 < Q 2 < 14 GeV 2 , 14 GeV 2 < Q 2 < 20 GeV 2 and 20 GeV 2 < 
Q 2 < 34 GeV 2 . Each data set is subdivided into bins of x v ; s of similar statistics, as listed in 
Table 2. A Bayesian unfolding procedure [20] is used to relate the measured x v i S to the true 
value of x and to correct the data for the detector acceptance and efficiency. This procedure 
is applied using, in turn, the PYTHIA and TWOGAM generators. The average of the cross 
sections obtained in the two cases is retained. The correlation between the generated value of 
x and x vis is similar for the two models. The one obtained with PYTHIA is shown in Figure 
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4b. The cross sections measured for each x interval of average value (x) and for each value 
of y/s are given in Table 2 with their statistical and systematic uncertainties. The bin-to-bin 
correlation matrices obtained with PYTHIA and TWOGAM for y/s = 189 GeV are shown in 
Tables 3 and 4, respectively. Similar matrices are obtained for the other values of y/s. 



4.2 Systematic Uncertainties 

The systematic uncertainties on the cross sections are estimated for each data sample, for each 
x bin and for each Q 2 interval. Three main sources of systematic uncertainties are considered: 
the selection procedure, the trigger efficiency and the Monte Carlo model. Their average effects 
over the full data sample are listed in Table 5. 

The uncertainties from the selection procedure are estimated by varying the selection cuts. 
The fiducial value of the polar angle in the low polar-angle calorimeter is varied from 0.0325 rad 
to 0.0360 rad and from 0.0637 rad to 0.060 rad. These changes result in a Q 2 -dependent 
uncertainty, as the highest and lowest values of the Q 2 are affected by the fiducial cut. The cut 
on E tag is varied from 0.70-Ebeam to 0. 65£beam and 0.75-Ebeam- The anti-tag cut is changed from 
0.20i?bcam to 0.15i?bcam and 0.25£beam- The numbers of particles is varied from four to three 
and five. The cut on the total energy in the calorimeters relative to \/s is varied from 0.40 to 
0.35 and 0.45. The uncertainty on the trigger efficiency, as determined from the data, is 1.5%. 
An additional uncertainty comes from the limited Monte Carlo statistics. 

The discrepancies of the results obtained with the PYTHIA and TWOGAM generators are 
considered as systematic uncertainties related to the Monte Carlo modelling. This difference 
is due to the calculated acceptance as well as to the unfolding procedure. An additional 
contribution to this modelling uncertainty is evaluated by repeating the analysis doubling or 
halving the VDM component of the TWOGAM generator and is found to be negligible. 
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4.3 Extraction of F 2 

The photon structure function /a is derived from the cross section of Equation 1 and the 
target-photon flux calculated by the program GALUGA [21] as: 



F2(x,Q 2 )/a = 



Acr meas (e + e — > e + e hadrons) 



Ao"Gai U ga(e + e — > e+e hadrons) 



The program calculates the theoretical value Aacaiuga in the given Q 2 and x range setting 
F 7 = 1 and F2 to the QPM value [22]. In practice the F2 contribution to the cross section 
is smaller than 1%, due to the small value of y. The running of the fine-structure constant 
with Q 2 is included. A GVDM form factor is used in the calculation for the target photon 
virtuality whose average value is of the order of 0.07 GeV 2 . The low polar-angle calorimeter 
acceptance for the tagged and the anti-tagged electron and the W 77 > 4 GeV requirement are 
taken into account. The uncertainty on Aacaiuga, as estimated by comparing the GVDM to a 
p form factor, is 2%. 

The contribution of radiative corrections to the cross section is evaluated by using the 
program RADCOR [23], which includes initial and final state radiation for the reaction e + e~ — > 
e + e~fi + fi~. The corrections are mainly due to initial state radiation from the electron scattered 
at large angle. Final state radiation is detected together with the scattered electron due to the 
finite granularity of the calorimeter. Initial state radiation from the electron producing the 
quasi-real target photon is very small. The calculations are performed at the generator level 
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using the Q 2 from the electron variables and W 77 from the muon pair. The measured F^jcx. is 
multiplied by the ratio, 71, of the non-radiative and the total cross section, shown in Tables 6 
and 7 for different values of x and Q 2 . 

The F2 values, averaged over the x intervals, are first obtained for each individual data 
set. The results are statistically compatible and, consequently, a weighted average of is 
calculated for the Q 2 ranges with average values 12.4 GeV 2 , 16.7 GeV 2 and 25.5 GeV 2 . This 
procedure is applied to data unfolded separately with PYTHIA and TWOGAM and the two 
different values are shown in Figure 5. Their average value is given in Table 6 and in Figure 6 
for each x interval of expected average value (x). In addition to the systematic uncertainty on 
the cross section, presented in Table 5, two systematic uncertainties are further considered in 
the extraction of F 2 7 . The first uncertainty, of 2%, is associated to the GALUGA calculations. 
The second uncertainty, also of 2%, covers the estimation of the initial-state and final-state 
radiative corrections. The uncertainty on initial-state radiation is assessed by changing the 
angular and momentum criteria which separate soft from hard photons in the Monte Carlo 
programs. The uncertainty on final-state radiation is estimated by varying the cone angle of 
the calorimeter for which final state radiation is detected together with the scattered electron. 

A comparison of the data with the existing parameterisations as obtained with the PDFLIB 
library [24] shows that our data are not well described by the leading-order parton density 
functions. In Figure 6 the data are compared with the predictions of the leading- and higher- 
order parton density functions GRV [25] and the higher-order parton density functions CJK [10]. 
The best agreement is found for the higher-order GRV [25] predictions. In all cases four quarks, 
u, d, s and c are used. The pure QPM prediction for 77 — > qq is also indicated, as calculated 
by using GALUGA with a mass of 0.32 GeV for the u and d quarks, 0.5 GeV for the s quark 
and 1.4 GeV for the c quark. It is clearly insufficient to describe the data. 

4.4 Q 2 -evolution of F£ 

The (^-evolution of F% , is studied in four x bins, x = 0.01 — 0.1, x = 0.1 — 0.2, x = 0.2 — 0.3 
and x = 0.3 — 0.5 and the results are given in Table 7. In Figure 7 the F^/ot values are 
presented for the lowest x bin and for a combined bin x — 0.1 — 0.5, together with our previous 
results [6]. Corrections for radiative effects are applied. The new measurements at 11 GeV 2 < 
Q 2 < 14 GeV 2 , 14 GeV 2 < Q 2 < 20 GeV 2 and 20 GeV 2 < Q 2 < 34 GeV 2 are in good agreement 
with our previous results. The expected linear growth with In Q 2 is observed in both x intervals. 
The function a + b\n(Q 2 / GeV 2 ) is fitted to the data, taking into account the total uncertainty 
calculated from the quadratic sum of statistical and systematic uncertainties. The fit results 
are: a = 0.141 ± 0.007 and b = 0.060 ± 0.005 for x = 0.01 - 0.1 with a confidence level of 44% 
and a = 0.05 ± 0.11 and b = 0.13 ± 0.04 for x = 0.1 - 0.5 with a confidence level of 71%. 

The predictions of the leading- and higher-order parton density functions GRV [25] and the 
higher-order parton density functions CJK [10] are also indicated in Figure 7. The evolution 
is different for the different models and the data are better described by the higher-order GRV 
model. 

5 Conclusions 

The photon structure function F^ is measured at LEP with the L3 detector at centre-of-mass 
energies 189 < y/s < 209 GeV in the Q 2 range 11 GeV 2 < Q 2 < 34 GeV 2 and the x range 
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0.006 < x < 0.556. The data are better reproduced by the higher-order parton density function 
of GRV than by other parton distribution functions determined from the low energy data. 

Combining the present results with previous L3 measurements, the Q 2 evolution is studied 
from 1.5 GeV 2 to 120 GeV 2 in the \ow-x region, 0.01 < x < 0.1, and from 12.4 GeV 2 to 
225 GeV 2 in the higher-x region, 0.1 < x < 0.5. The measurements at different centre-of-mass 
energies are consistent and the In Q 2 evolution of F 2 7 is clearly confirmed. 
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<Vi)(GeV) 


Apb" 1 ) 


Events 


e 


t-e- - 


-> e + e t + t 


e+e -> qq(7) 


Purity (%) 


Q 


) 2 , 


= 11 


- 14 GeV 2 






189 


171.8 


1884 








107.6 


5.7 


94 


194 


111.4 


1197 








76.3 


3.3 


93 


200 


109.3 


1187 








74.9 


3.7 


93 


206 


215.6 


2418 








129.6 


7.7 


94 








Q 2 - 


= 14 


- 20 GeV 2 






189 


171.8 


2046 








128.6 


9.7 


93 


194 


111.4 


1347 








91.2 


5.4 


93 


200 


109.3 


1359 








89.5 


4.9 


93 


206 


215.6 


2886 








177.4 


8.7 


94 








Q 2 - 


= 20 


- 34 GeV 2 






189 


171.8 


1922 








143.9 


8.1 


92 


194 


111.4 


1331 








103.6 


5.8 


92 


200 


109.3 


1287 








101.6 


6.7 


92 


206 


215.6 


2859 








202.5 


12.7 


92 



Table 1: The average e + e centre-of-mass energies, (y/s), and the corresponding luminosities for 
the four data samples together with the number of selected events in the Q 2 intervals 1 1 GeV 2 < 
Q 2 < 14 GeV 2 , 14 GeV 2 < Q 2 < 20 GeV 2 and 20 GeV 2 < Q 2 < 34 GeV 2 . The numbers of 
expected background events from the e + e~ — > e + e~r + r~ and e + e~ — > qq(7) processes and the 
signal purity are also listed. 
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189 GeV 


194 GeV 


200 GeV 


zUo Lrev 












Q 2 = n- 


- 14 GeV^ 










Ax 


range 


/™\ 
(x) 


Aa/Ax (pb) 


Aa/Ax (pb) 


Aa/Ax (pb) 


A /A / 1 ~\ 

Aa/Ax (pb) 


0.006 


— 0.023 


0.013 


103.6 ± 8.5 ± 14.5 


108.2 ±10.8 ±15.4 


106.4 ± 10.2 ± 10.4 


1 1 r o 1 TO 1 1 /"* 4 

115. o ± 7.o ± lb. 4 


0.023 


- 0.040 


0.031 


63.4 ± 5.3 ± 


4.3 


67.0 ± 


6.5 ± 8.4 


63.9 ± 


6.3 ± 


5.2 


an c: I A c I T T 
by.b ± 4.0 ± 7.7 


0.040 


- 0.060 


0.050 


52.0 ± 4.0 ± 


2.9 


55.7 ± 


5.1 ± 7.8 


53.4 ± 


5.5 ± 


5.3 


r n n i o c i r a 

52.0 ± 3.6 ± 5.4 


0.060 


- 0.090 


0.075 


45.0 ± 3.4 ± 


2.2 


47.7 ± 


4.4 ± 6.1 


44.0 ± 


4.0 ± 


4.4 


/IOO 1 O O 1 f o 

43.8 ± 2.8 ± 5.3 


0.090 


- 0.120 


0.10 


40.2 ± 2.9 ± 


1.8 


40.8 ± 


3.9 ± 5.3 


39.3 ± 


3.6 ± 


4.3 


on /? 1 O o 1 C 1 

39.6 ± 2.8 ± 5.1 


0.120 


— 0.160 


0.14 


37.9 ± 2.8 ± 


2.3 


37.4 ± 


3.4 ± 5.3 


36.9 ± 


3.4 ± 


4.2 


o t n i n r i a a 

37.9 ± 2.5 ± 4.4 


0.160 


- 0.205 


0.18 


34.8 ± 2.5 ± 


1.3 


33.9 ± 


2.9 ± 4.7 


33.0 ± 


2.9 ± 


4.3 


o r O 1 O A 1 A O 

35.8 ± 2.4 ± 4.3 


0.205 


- 0.260 


0.23 


33.3 ± 2.4 ± 


1.3 


31.9 ± 


3.1 ± 4.8 


31.6 ± 


2.8 ± 


4.0 


oo 1 C\ C\ 1 A A 

32.6 ± 2.0 ± 4.4 


0.260 


- 0.330 


0.29 


29.8 ± 2.1 ± 


1.7 


29.1 ± 


2.9 ± 4.0 


28.3 ± 


2.7 ± 


6. 7 


o r\ T i -i r\ i a A 

30.7 ± 1.9 ± 4.4 


0.330 


- 0.400 


0.36 


29.9 ± 2.1 ± 


1.3 


27.5 ± 


2.7 ± 4.0 


25.5 ± 


2.6 ± 


Q O 

o.o 


on r i i a i on 

29.5 ± 1.8 ± 3.9 












Q 2 = 14 - 


- 20 GeV^ 










0.006 


— 0.023 


0.013 


97.0 ± 8.5 ± 12.7 


100.4 ±10.2 ±9.8 


102.4 ± 10.4 ± 


11.1 


113.0 ± 7.7 ± 11.5 


0.023 


n n ac\ 

— 0.040 


n Hoi 

0.031 


59.3 ± 4.9 ± 


3.6 


64.0 ± 


6.5 ±4.9 


64.2 ± 


7.1 ± 


6.2 


170 1 tr 1 To 

72.2 ± 0.2 ± 7.0 


0.040 


- 0.060 


0.050 


49.6 ± 4.3 ± 


5.7 


50.8 ± 


5.2 ±5.0 


53.4 ± 


5.3 ± 


4.3 


rr O I O a I T 

00.0 ± 0.0 ± 2. 7 


0.060 


- 0.090 


0.075 


40.8 ± 2.9 ± 


1.7 


42.5 ± 


3.9 ±1.9 


44.5 ± 


4.2 ± 


2.6 


at 1 on 1 01 

47.8 ± 2.9 ± 3.1 


0.090 


- 0.120 


0.10 


36.1 ± 2.9 ± 


1.9 


38.0 ± 


3.6 ±2.3 


37.9 ± 


3.6 ± 


3.2 


a ~t r* 1 A 1 on 

41.5 ± 2.4 ± 2.0 


0.120 


- 0.160 


0.14 


32.5 ± 2.2 ± 


1.5 


35.7 ± 


3.0 ± 1.6 


36.7 ± 


3.3 ± 


2.1 


O T T 1 O O 1 1 1-7 

37.7 ± 2.3 ± 1.7 


0.160 


- 0.205 


0.18 


31.8 ± 2.3 ± 


1.3 


31.9 ± 


2.7 ± 1.7 


33.3 ± 


3.0 ± 


1.7 


O A T 1 O O 1 If 

34.7 ± 2.2 ± 1.5 


0.205 


- 0.260 


0.23 


30.8 ± 2.1 ± 


1.0 


30.6 ± 


2.8 ±1.3 


32.2 ± 


2.8 ± 


1.8 


A A 1 n 1 1 r 

34.4 ± 2.0 ± 1.5 


0.260 


- 0.330 


0.29 


27.7 ± 2.0 ± 


1.0 


27.6 ± 


2.3 ±1.2 


30.0 ± 


2.3 ± 


1.7 


0101 "i n 1 10 

31.2 ± 1.9 ± 1.2 


0.330 


- 0.400 


0.36 


28.1 ± 2.1 ± 


1.1 


26.8 ± 


2.4 ±1.3 


29.4 ± 


2.8 ± 


1.3 


on n 1 1 n 1 1 n 

29.0 ± 1.9 ± 1.0 


0.400 


- 0.467 


0.42 


26.7 ± 2.1 ± 


2.1 


25.2 ± 


2.5 ± 1.8 


26.8 ± 


2.7 ± 


1.7 


T 1 1 n 1 11 

27.8 ± 1.9 ± 1.1 












Q 2 = 20 - 


- 34 GeV^ 










r\ noo 

0.023 


n n a n 

— 0.040 


n no i 

0.031 


48.6 ±5.4 ±3.6 


54.3 ±6.5 ±3.3 


37.4 ±3.6 ±2.4 


39.5 ± 2.6 ± 2.0 


0.040 


- 0.060 


0.050 


41.9 ±4.1 ±4.1 


36.0 ±3.4 ±2.4 


37.4 ±3.6 ±2.4 


on r 1 1 n 
oy.o ± 2.0 ± Z.U 


0.060 


- 0.090 


0.075 


35.8 ±2.8 ±1.7 


36.0 ±3.4 ±2.4 


37.4 ±3.6 ±2.4 


39.5 ±2.6 ±2.0 


0.090 


-0.120 


0.10 


32.8 ±2.7 ±2.2 


32.1 ±3.3 ± 1.4 


32.4 ± 3.3 ± 1.9 


36.3 ±2.6 ±1.3 


0.120 


-0.160 


0.14 


29.1 ±2.5 ± 1.7 


30.6 ±2.8 ±2.4 


32.1 ±3.0 ±2.5 


33.6 ±2.2 ±1.2 


0.160 


- 0.205 


0.18 


27.7 ±2.1 ±1.1 


28.1 ±2.6 ± 1.1 


28.6 ±2.6 ±1.7 


32.2 ±2.0 ±1.2 


0.205 


- 0.260 


0.23 


25.0 ±1.7 ±1.0 


24.2 ±2.0 ±1.0 


25.0 ± 2.2 ± 1.3 


27.4 ±1.6 ±1.2 


0.260 


- 0.330 


0.29 


25.0 ±1.7 ±1.0 


24.2 ±2.0 ±1.0 


25.0 ± 2.2 ± 1.3 


27.4 ± 1.6 ± 1.2 


0.330 


- 0.400 


0.36 


24.5 ± 1.8 ± 1.2 


23.3 ± 2.2 ± 1.0 


24.3 ± 2.2 ± 1.4 


25.6 ± 1.6 ±0.9 


0.400 


- 0.467 


0.42 


23.6 ± 1.8 ± 1.5 


22.4 ± 2.0 ± 1.1 


23.5 ±2.4 ±1.2 


24.5 ± 1.7 ± 1.0 


0.467 


- 0.556 


0.49 


22.8 ±1.6 ±2.5 


22.4 ±2.1 ±2.9 


24.0 ±2.3 ±2.3 


24.5 ± 1.5 ± 1.5 



Table 2: Cross sections Aa/Ax as a function of x for the reaction e + e~— > e + e~hadrons for the 
four average values of y/s, in three Q 2 intervals. The first uncertainty is statistical, the second 
systematic. 
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2 - 1 1 


- 14 GcV^ 










u.uuu 


093 
u.uzo 


OzLO 
U.U^iU 


060 
u.uuu 


U.Uju 


n i on 


n 1 (\a 

U. 1DU 


o 90 s ; 


960 
u.zuu 


330 


























023 


040 


060 


090 


1 90 


IfiO 

U. 1UU 


\J.Zi\JO 


960 


330 


400 


006 — 023 


1.00 




















023 - 040 


0.92 


1.00 


















0.040- 0.060 


0.75 


0.92 


1.00 
















0.060- 0.090 


0.55 


0.79 


0.96 


1.00 














0.090- 0.120 


0.45 


0.67 


0.85 


0.94 


1.00 












0.120-0.160 


0.39 


0.62 


0.78 


0.89 


0.98 


1.00 










0.160- 0.205 


0.28 


0.50 


0.67 


0.80 


0.95 


0.97 


1.00 








0.205 - 0.260 


0.25 


0.45 


0.59 


0.69 


0.84 


0.87 


0.94 


1.00 






0.260- 0.330 


0.22 


0.38 


0.49 


0.59 


0.79 


0.82 


0.91 


0.99 


1.00 




0.330- 0.400 


0.17 


0.36 


0.50 


0.60 


0.77 


0.78 


0.86 


0.98 


0.97 


1.00 













Q 2 = 


14 - 20 GeV^ 












0.006 


0.023 


0.040 


0.060 


0.090 


0.120 


0.160 


0.205 


0.260 


0.330 


0.400 


x range 


























0.023 


0.040 


0.060 


0.090 


0.120 


0.160 


0.205 


0.260 


0.330 


0.400 


0.467 


0.006 - 0.023 


1.00 






















0.023 - 0.040 


0.87 


1.00 




















0.040 - 0.060 


0.74 


0.95 


1.00 


















0.060 - 0.090 


0.55 


0.82 


0.93 


1.00 
















0.090 - 0.120 


0.54 


0.79 


0.89 


0.98 


1.00 














0.120-0.160 


0.36 


0.63 


0.73 


0.89 


0.93 


1.00 












0.160 - 0.205 


0.31 


0.60 


0.71 


0.85 


0.89 


0.98 


1.00 










0.205 - 0.260 


0.24 


0.49 


0.59 


0.75 


0.81 


0.95 


0.95 


1.00 








0.260 - 0.330 


0.18 


0.40 


0.50 


0.64 


0.71 


0.87 


0.88 


0.97 


1.00 






0.330 - 0.400 


0.21 


0.40 


0.47 


0.56 


0.65 


0.79 


0.79 


0.91 


0.97 


1.00 




0.400 - 0.467 


0.11 


0.31 


0.37 


0.47 


0.55 


0.71 


0.71 


0.84 


0.91 


0.96 


1.00 













Q 2 = 


20 - 34 GcV A 












0.023 


0.040 


0.060 


0.090 


0.120 


0.160 


0.205 


0.260 


0.330 


0.400 


0.467 


x range 


























0.040 


0.060 


0.090 


0.120 


0.160 


0.205 


0.260 


0.330 


0.400 


0.467 


0.556 


0.023 - 0.040 


1.00 






















0.040 - 0.060 


0.92 


1.00 




















0.060 - 0.090 


0.87 


0.96 


1.00 


















0.090 - 0.120 


0.80 


0.88 


0.97 


1.00 
















0.120-0.160 


0.66 


0.75 


0.89 


0.96 


1.00 














0.160 - 0.205 


0.62 


0.71 


0.82 


0.91 


0.97 


1.00 












0.205 - 0.260 


0.58 


0.64 


0.76 


0.86 


0.93 


0.98 


1.00 










0.260 - 0.330 


0.47 


0.54 


0.66 


0.77 


0.88 


0.95 


0.97 


1.00 








0.330 - 0.400 


0.43 


0.49 


0.61 


0.71 


0.81 


0.88 


0.92 


0.97 


1.00 






0.400 - 0.467 


0.38 


0.45 


0.55 


0.63 


0.73 


0.78 


0.83 


0.91 


0.98 


1.00 




0.467-0.556 


0.39 


0.44 


0.54 


0.63 


0.73 


0.78 


0.84 


0.91 


0.98 


0.99 


1.00 



Table 3: Correlation matrices obtained with the PYTHIA Monte Carlo for the data at y/s = 
189 GeV for the three Q 2 intervals. 
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n 


2 = 11- 


- 14 GcV z 










0.006 


0.023 


0.040 


0.060 


0.090 


0.120 


0.160 


0.205 


0.260 


0.330 




0.023 


0.040 


0.060 


0.090 


0.120 


0.160 


0.205 


0.260 


0.330 


0.400 


0.006-0.023 


1.00 




















0.023-0.040 


0.91 


1.00 


















0.040-0.060 


0.79 


0.96 


1.00 
















0.060-0.090 


0.63 


0.86 


0.95 


1.00 














0.090-0.120 


0.52 


0.74 


0.86 


0.95 


1.00 












0.120-0.160 


0.46 


0.66 


0.78 


0.90 


0.98 


1.00 










0.160-0.205 


0.36 


0.57 


0.69 


0.82 


0.92 


0.97 


1.00 








0.205-0.260 


0.27 


0.45 


0.56 


0.69 


0.80 


0.88 


0.97 


1.00 






0.260-0.330 


0.24 


0.42 


0.52 


0.65 


0.73 


0.83 


0.93 


0.98 


1.00 




0.330-0.400 


0.19 


0.36 


0.45 


0.58 


0.64 


0.74 


0.86 


0.93 


0.98 


1.00 













o 2 - 


14 - 20 GcV^ 












0.006 


0.023 


0.040 


0.060 


0.090 


0.120 


0.160 


0.205 


0.260 


0.330 


0.400 




0.023 


0.040 


0.060 


0.090 


0.120 


0.160 


0.205 


0.260 


0.330 


0.400 


0.467 


0.006-0.023 


1.00 






















0.023-0.040 


0.90 


1.00 




















0.040-0.060 


0.76 


0.95 


1.00 


















0.060-0.090 


0.62 


0.86 


0.96 


1.00 
















0.090-0.120 


0.50 


0.74 


0.87 


0.96 


1.00 














0.120-0.160 


0.43 


0.66 


0.78 


0.90 


0.97 


1.00 












0.160-0.205 


0.36 


0.56 


0.68 


0.81 


0.90 


0.96 


1.00 










0.205-0.260 


0.29 


0.45 


0.57 


0.70 


0.78 


0.87 


0.96 


1.00 








0.260-0.330 


0.23 


0.38 


0.49 


0.60 


0.68 


0.77 


0.90 


0.97 


1.00 






0.330-0.400 


0.17 


0.32 


0.43 


0.56 


0.63 


0.73 


0.85 


0.95 


0.98 


1.00 




0.400-0.467 


0.16 


0.28 


0.38 


0.48 


0.55 


0.64 


0.77 


0.88 


0.95 


0.98 


1.00 




































Q 2 = 


20 - 34 GeV z 












0.023 


0.040 


0.060 


0.090 


0.120 


0.160 


0.205 


0.260 


0.330 


0.400 


0.467 


x range 


0.040 


0.060 


0.090 


0.120 


0.160 


0.205 


0.260 


0.330 


0.400 


0.467 


0.556 


0.023-0.040 


1.00 






















0.040-0.060 


0.92 


1.00 




















0.060-0.090 


0.85 


0.95 


1.00 


















0.090-0.120 


0.73 


0.87 


0.96 


1.00 
















0.120-0.160 


0.62 


0.77 


0.88 


0.96 


1.00 














0.160-0.205 


0.52 


0.66 


0.78 


0.88 


0.96 


1.00 












0.205-0.260 


0.45 


0.61 


0.71 


0.81 


0.91 


0.98 


1.00 










0.260-0.330 


0.35 


0.50 


0.59 


0.69 


0.79 


0.90 


0.96 


1.00 








0.330-0.400 


0.34 


0.45 


0.54 


0.62 


0.69 


0.78 


0.85 


0.93 


1.00 






0.400-0.467 


0.25 


0.37 


0.46 


0.56 


0.66 


0.78 


0.87 


0.96 


0.98 


1.00 




0.467-0.556 


0.25 


0.36 


0.44 


0.52 


0.61 


0.72 


0.81 


0.90 


0.97 


0.98 


1.00 



Table 4: Correlation matrices obtained with the TWOGAM Monte Carlo for the data at 
y/s = 189 GeV for the three Q 2 intervals. 
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Source of systematic uncertainties 


Uncertainty in % 


Tagging calorimeter polar angle 


0.9-8.0 


Tagging calorimeter energy 


1.6-2.5 


Anti-tag energy 


0.4 


Number of particles 


0.2-2.6 


Total energy in the calorimeters 


0.2 


Trigger efficiency 


1.5 


Monte Carlo statistics 


0.2-0.5 


Model dependence 


0.7-9.9 



Table 5: Systematic uncertainties on the measured cross sections. 
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x range 


(x) 


K 


FZ/a 






Q 2 = ll 


- 14 GeV^ 


0.006 - 


0.023 


0.013 


0.92 


0.302 


±0.013 ±0.026 ± 0.029 


0.023 - 


0.040 


0.031 


0.90 


0.245 


±0.011 ±0.021 ± 0.010 


0.040 - 


0.060 


0.050 


0.88 


0.257 


±0.011 ±0.023 ± 0.012 


0.060 - 


0.090 


0.075 


0.90 


0.296 


±0.012 ±0.028 ± 0.009 


0.090 - 


0.120 


0.10 


0.89 


0.315 


±0.013 ±0.032 ± 0.007 


0.120- 


0.160 


0.14 


0.90 


0.365 


±0.015 ±0.038 ± 0.008 


0.160- 


0.205 


0.18 


0.88 


0.399 


±0.017 ±0.043 ± 0.007 


0.205 - 


0.260 


0.23 


0.89 


0.441 


±0.018 ±0.049 ± 0.011 


0.260 - 


0.330 


0.29 


0.88 


0.483 


±0.020 ±0.054 ± 0.013 


0.330 - 


0.400 


0.36 


0.89 


0.536 


±0.023 ±0.059 ± 0.025 






Q 2 = 14 


- 20 GeV 2 


0.006 - 


0.023 


0.013 


0.93 


0.310 


±0.014 ±0.020 ± 0.028 


0.023 - 


0.040 


0.031 


0.88 


0.258 


±0.012 ±0.014 ± 0.018 


0.040 - 


0.060 


0.050 


0.90 


0.275 


±0.012 ±0.014 ± 0.017 


0.060 - 


0.090 


0.075 


0.90 


0.288 


±0.011 ±0.014 ± 0.008 


0.090 - 


0.120 


0.10 


0.89 


0.316 


±0.013 ±0.015 ± 0.014 


0.120 - 


0.160 


0.14 


0.90 


0.337 


±0.013 ±0.016 ± 0.007 


0.160- 


0.205 


0.18 


0.89 


0.381 


±0.015 ±0.018 ± 0.006 


0.205 - 


0.260 


0.23 


0.88 


0.424 


±0.017 ±0.019 ± 0.008 


0.260 - 


0.330 


0.29 


0.88 


0.471 


±0.018 ±0.020 ± 0.009 


0.330 - 


0.400 


0.36 


0.87 


0.510 


±0.021 ±0.022 ± 0.007 


0.400 - 


0.467 


0.42 


0.88 


0.551 


±0.024 ±0.024 ± 0.026 






Q 2 = 20 


- 34 GeV 2 


0.023 - 


0.040 


0.031 


0.89 


0.317 


±0.017 ± 0.016 ±0.017 


0.040 - 


0.060 


0.050 


0.89 


0.293 


±0.015 ± 0.014 ±0.010 


0.060 - 


0.090 


0.075 


0.89 


0.314 


±0.013 ± 0.015 ±0.012 


0.090 - 


0.120 


0.10 


0.88 


0.338 


±0.016 ± 0.016 ±0.018 


0.120 - 


0.160 


0.14 


0.88 


0.384 


±0.017 ± 0.018 ±0.017 


0.160- 


0.205 


0.18 


0.88 


0.404 


±0.017 ± 0.018 ±0.004 


0.205 - 


0.260 


0.23 


0.88 


0.446 


±0.018 ± 0.020 ±0.009 


0.260 - 


0.330 


0.29 


0.87 


0.488 


±0.019 ± 0.022 ±0.006 


0.330 - 


0.400 


0.36 


0.87 


0.557 


±0.023 ± 0.025 ±0.010 


0.400 - 


0.467 


0.42 


0.87 


0.611 


±0.027 ± 0.029 ±0.015 


0.467- 


0.556 


0.49 


0.87 


0.683 


±0.028 ± 0.030 ±0.062 



Table 6: Measured values of F^ jot and the applied radiative correction factors, 71, in bins of x 
for the three Q 2 ranges. The first uncertainty is statistic, the second systematic and the third 
is due to model dependence. 
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Q 2 range 


(Q 2 ) 


K 


Fq/ot 






x = 


-- 0.01-0.1 


11 - 14 


12.4 


0.89 


0.278 ± 0.006 ± 0.028 ± 0.013 


14 - 20 


16.7 


0.89 


0.287 ± 0.006 ± 0.015 ± 0.015 


20 - 34 


25.5 


0.88 


0.316 ± 0.008 ± 0.016 ± 0.013 






x - 


= 0.1-0.2 


11 - 14 


12.4 


0.88 


0.377 ± 0.010 ± 0.039 ± 0.008 


14-20 


16.7 


0.88 


0.355 ± 0.009 ± 0.017 ± 0.005 


20-34 


25.5 


0.88 


0.399 ± 0.011 ± 0.019 ± 0.010 






x - 


= 0.2-0.3 


11 - 14 


12.4 


0.88 


0.464 ± 0.015 ± 0.051 ± 0.009 


14 - 20 


16.7 


0.88 


0.442 ± 0.013 ± 0.020 ± 0.003 


20 - 34 


25.5 


0.87 


0.477 ± 0.015 ± 0.023 ± 0.013 






x - 


= 0.3-0.5 


11-14 


12.4 


0.89 


0.544 ± 0.017 ± 0.061 ± 0.019 


14-20 


16.7 


0.87 


0.545 ± 0.014 ± 0.024 ± 0.012 


20-34 


25.5 


0.87 


0.594 ± 0.015 ± 0.029 ± 0.022 



Table 7: The values of F^jct in bins of Q 2 for four x ranges together with the radiative correction 
factor. The first uncertainty is statistic, the second systematic and the third is due to model 
dependence. 
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Figure 1: Distribution of the highest energy clusters in the forward electromagnetic calorimeters 
for a) the tagged electron side and b) for the opposite side, c) Total energy in the central 
calorimeters, d) The visible mass of the hadronic final state. All distributions are presented 
after all other cuts are applied. The backgrounds from annihilation and two-photon leptonic 
events are indicated as shaded areas and added to the expectations of the PYTHIA, PHOJET 
and TWOGAM generators. The arrows indicate the position of the cuts. 
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Figure 2: Q 2 distribution of the selected events for the four average ^Js ranges. 
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Figure 3: Distribution of the visible mass of the two-photon system and of the x v - ls for all 
selected events compared with Monte Carlo predictions. The backgrounds from annihilation 
and two-photon leptonic events are indicated as shaded areas and added to the expectations of 
the PYTHIA, PHOJET and TWOGAM generators. 
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Figure 4: a) The detector acceptance and selection efficiency, e, obtained by the PYTHIA and TWOGAM generators, b) Comparison 
of the reconstructed and generated value of x for the PYTHIA Monte Carlo at y/s = 189 GeV for different values of Q 2 . The mean 
observed value and the standard deviation of x v i S are plotted for events generated in a given x bin. 
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Figure 5: The photon structure function F^/a obtained with PYTHIA and TWOGAM. Only 
the statistical uncertainties are shown. For clarity, the symbols corresponding to the two Monte 
Carlo generators are slightly offset. The measurements are correlated as indicated in Tables 3 
and 4 
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Figure 6: The photon structure function F2 /a as a function of x for the three Q 2 intervals, 
with statistical and systematic uncertainties. The former are indicated by the inner error bars. 
The new data are presented together with the previous results at y/s = 183 GeV [6]. The 
predictions of the higher-order parton density functions GRV and CJK are shown as well as 
the leading-order predictions of the GRV. The changes in slope of the CJK predictions are due 
to the c-quark threshold. The QPM predictions for 77 — > qq are also shown. 
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Figure 7: Evolution of the photon structure function F^/a as a function of Q 2 for two x 
intervals. The results of a fit to the data of the function a + b(\nQ 2 / GeV 2 ) are shown together 
with the predictions of the higher-order parton density functions GRV and CJK as well as the 
leading-order predictions of GRV. 
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